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Abstract 

We discuss the production of particles in relativistic heavy-ion collisions 
through the mechanism of massive bremsstrahlung, in which massive mesons 
are emitted during rapid nucleon acceleration. This mechanism is described 
within the framework of classical hadrodynamics for extended nucleons, corre- 
sponding to nucleons of finite size interacting with massive meson fields. This 
new theory provides a natural covariant microscopic approach to relativistic 
heavy-ion collisions that includes automatically spacetime nonlocality and retar- 
dation, nonequilibrium phenomena, interactions among all nucleons, and par- 
ticle production. Inclusion of the finite nucleon size cures the difficulties with 
preacceleration and runaway solutions that have plagued the classical theory of 
self- inter acting point particles. For the soft reactions that dominate nucleon- 
nucleon collisions, a significant fraction of the incident center-of-mass energy is 
radiated through massive bremsstrahlung. In the present version of the theory, 
this radiated energy is in the form of neutral scalar (a) and neutral vector (u) 
mesons, which subsequently decay primarily into pious with some photons also. 
Additional meson fields that are known to be important from nucleon-nucleon 
scattering experiments should be incorporated in the future, in which case the 
radiated energy would also contain isovector pseudoscalar (vr"*", vr", vr*^), isovec- 
tor scalar {6~^ , 6~ , 6^), isovector vector p~ , p^), and neutral pseudoscalar 
(rj) mesons. 



1. Introduction 

Many particles are produced in a typical relativistic heavy-ion collision.S'i For 
the production of these particles, we would like to discuss the mechanism of massive 
bremsstrahlung, in which massive mesons are emitted during rapid nucleon accelera- 
tion. This mechanism is described within the framework of classical hadrodynamics 
for extended nucleons, corresponding to nucleons of finite size interacting with massive 
meson fields. This approach, which satisfies a priori the physical conditions that exist 



Fig. 1: Slice through the center of a nucleon. The circle indicates the location of the root-mean-square 
radius, where the exponentially decreasing mass density is only 3% of its central value. 



at relativistic energies, is manifestly Lorentz-covariant and allows for nonequilibrium 
phenomena, interactions among all nucleons, and particle production. 

Although the nucleon is a composite particle made up of three valence quarks plus 
additional sea quarks and gluons, when nucleons collide at very high energies, only a 
few rare events correspond to the head-on or hard collisions between the individual 
quarks and/or gluons. Whereas the underlying quark-gluon structure of the nucleon 
is of crucial importance for describing particle production in hard collisions, such 
collisions are nevertheless extremely rare, typically one in a billion.! The vast majority 
of events correspond to soft collisions not involving individual quarks or gluons. For 
describing particle production in such events, an appealing idea is to regard the nucleon 
as a single extended object interacting with other nucleons through the conventional 
exchange of mesons (whose underlying quark-antiquark composition is ignored). 



the proton charge density is approximately exponential in shape,u with a root-mean- 
square radius of 0.862 ± 0.012 fm. Although many questions remain concerning the 
relationship between the proton charge density and the nucleon mass density,! it should 
be a fairly accurate approximation to regard them as equal. We therefore take the 
nucleon mass density to be 




s involving elastic electron scattering off protons have determined that 
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with /i = ^/T2/-Rrms and -Rrms = 0.862 fm. We show in fig. |1| a gray-scale plot of the 
mass density through the center of a nucleon calculated according to this exponential, 
with the root-mean-square radius indicated by a circle. 

The physical input underlying our new approach consists of Lorentz invariance 
(which includes energy and momentum conservation), nucleons of finite size interacting 
with massive meson fields, and the classical approximation applied in domains where 
it should be reasonably valid. At bombarding energies of many GeV per nucleon, 
the de Broglie wavelength of projectile nucleons is extremely small compared to all 
other length scales in the problem. In addition, the Compton wavelength of the 
nucleon is small compared to its radius, so that effects due to the intrinsic size of the 
nucleon dominate those due to quantum uncertainty. Finally, the angular momentum 
is typically several hundred and the radiated energy corresponds to several meson 
masses. The classical approximation for nucleon trajectories should therefore be valid, 
provided that the effect of the finite nucleon size on the equations of motion is taken 
into account. 

We describe in sect. 2 the present version of our theory, which includes the neutral 
scalar (a) and neutral vector (cu) meson fields. This permits a qualitative discussion 
of not only particle production through massive bremsstrahlung, but also such other 
physically relevant points as the effect of the finite nucleon size on the equations of 
motion and an inherent spacetime nonlocality that may be responsible for significant 
collective effects. The a and oj mesons that are produced will subsequently decay pri- 
marily into pions with some photons also. The resulting classical relativistic equations 
of motion are solved in sect. 3 for soft nucleon-nucleon collisions at piab = 14.6, 30, 60, 
100, and 200 GeV/c. Section 4 discusses the future incorporation of additional meson 
fields that are known to be important from nucleon-nucleon scattering experiments, 
including the isovector pseudoscalar (tt^, tt", 7r°), isovector scalar 5", 5°), isovec- 
tor vector (p"*", P", p°), and neutral pseudoscalar {rj). Further details are given in a 
series of papers,Q"Ej although not all of the equations appearing in some of the earlier 
publications are in their final form. 



2. Equations of motion 

Our action for extended, unexcited nucleons interacting with massive scalar and 
vector meson fields is 
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Vector field Interaction 

where Mq is the bare nucleon mass and qi is the four-position of the ith nucleon, 
whose trajectory is given by qi = qi{Ti). A dot represents the derivative with respect 
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to Tj. In the action the four- velocities are not constrained so that = 1 and Xj is 
not yet identified as the proper time; it is only in the equations of motion, which 
are derived as a result of the variation of /, that this is true. We use the metric 
g^'^ = diag(l, — 1, — 1, — 1), write four-vectors as = {q^,q) = {q^,q^,q^,q^), and 
use units in which h = c = 1. The scalar potential is denoted by 0, the four- vector 
potential by V, and the meson masses by ms,v The vector field strength tensor is 



the scalar source density is 
and the vector source density is 
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where p is the four-dimensional mass density of the nucleon, the spatial iiart of which 
we assume to be exponential in the nucleon's rest frame. The values@"lii that we 
have used here for the six physical constants appearing in our theory are nucleon mass 
M = 938.91897 MeV, scalar (a) meson mass rris = 550 MeV, vector [u) meson mass 
m-v = 781.95 MeV, scalar interaction strength gj^ = 7.29, vector interaction strength 
(yfv^ = 10.81, and nucleon r.m.s. radius Rrms = 0.862 fm. 

In ref. we have derived exact equations of motion for the above action in two 
limits: (1) relativistic point nucleons and (2) nonrelativistic extended nucleons. We 
then generalize covariantly to obtain relativistic equations of motion for extended 
nucleons, which can be written as 



M*n'^ = -u -u 4- 
The effective mass is given by 

AMsclf,i + fi-s^ext,* , 



M* = M* 
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The hadrostatic sclf-cnergies of the scalar and vector fields are denoted by Mg ^ and 
their logarithmic derivatives with respect to meson mass by Mg' ,^. The self- forces are 
given by 
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and the external forces by 
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In the above = 9^"" - s'j = - Qjirj - a), vj = gj(r,), v'j = qj{Tj - a), 

ttj = qj{Tj), kj = s'j ■ q'j, and the retarded proper time Tj is determined implicitly from 
the condition kj — 0, with a — 0. 

These equations are written in terms of the nucleon structure functions h and w 
and quantities derived from them. We define the interaction energy function to be 



(15) 



where R= |xi — X2I and r is the distance between the centers of the two particles. Here 
p(r) is the nucleon mass density normalized so that An r'^drp{r) = 1. A Laplace 
transform relates the structure function w and its derivative w'{a,r) = r~^dw/dr to 
the interaction energy W: 



w{a,r) = 2£(-^)[W^(m,r);2a] . 
The self-interaction structure function is 



h{a) = 3271^ da' (a'^ - a^) p{a + a')p{\a - a'\) , 



(16) 

(17) 
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Fig. 2: Calculated dependence of the center-of-mass scattering angle upon impact parameter for soft 
nucleon-nuclcon collisions at five incident laboratory momenta. 

with h'{a) = dh{a)/da, h{a) = Jq da'h{a'), and h{a) = h'{cr) — Qm^'^h{a). The 
first-order Bessel function of the first kind is denoted by Ji. 

These equations of motion, which are second-order, nonhnear, integrodifferential 
equations with four dimensions per particle, can be solved numerically without further 
approximation. In particular, we do not need to make either a mean-field approxima- 
tion, a perturbative expansion in coupling strength, or a superposition of two-body 
collisions. To solve them we use a fourth-order Adams-Moulton predictor-corrector 
algorithm with adaptive step sizes. The integrations over proper time are done with 
a special error-minimizing application of Lagrange's four-point (cubic) interpolation 
formulas. 

3. Radiated energy and other results for soft nucleon-nucleon collisions 

We now present some results obtained by solving our equations of motion for 
the soft collision of two nucleons at laboratory momentum piab = 14.6, 30, 60, 100, 
and 200 GeV/c. At three of these momenta substantial experimental data exist for 
heavy-ion collisions,^'! and at the remaining two momenta experimental data exist for 
proton-proton collisions. Ii3 We will concentrate our discussion here on such physically 
observable quantities as scattering angle, transverse momentum, and radiated energy 
in the center-of-mass system, in which frame the computations are performed. 

As shown in fig. ^, for a given incident momentum, the center-of-mass scattering 
angle for the dominating soft reactions described by our theory has a maximum value 
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Fig. 3: Calculated dependence of the transverse momentum upon impact parameter for soft nucleon- 
nucleon collisions at five incident laboratory momenta. 



ace at which the trans- 
^ has its maximum. At 



at a certain impact parameter and decreases to zero for both head-on and distant colh- 
sions. With increasing incident momentum in this range both the maximum angle and 
the impact parameter at which this maximum occurs decrease. For ultrarelativistic 
colhsions this impact parameter is approximately the dista: 
versely dominating static vector force for extended nucleonJH 
the other extreme of low incident momentum, the opposing scalar and vector forces are 
of similar magnitude and give rise for small impact parameter to the more complicated 
behavior of the double-dot-dashed curve in fig. ^. 

The transverse momentum has a related behavior, as shown in fig. ^. For a given 
incident momentum, the transverse momentum for soft reactions also has a maximum 
value at a certain impact parameter and decreases to zero for both head-on and dis- 
tant collisions. The maximum transverse momentum increases slowly with increasing 
incident momentum in this range, and the impact parameter at which this maximum 
occurs decreases. 

The center-of-mass radiated energy per nucleon for soft reactions shown in fig. ^ also 
has a maximum value at a certain impact parameter. However, this quantity decreases 
to a finite value for head-on collisions and to zero for distant collisions. The maximum 
center-of-mass radiated energy per nucleon increases strongly with increasing incident 
momentum. In the present version of the theory, this radiated energy will be in the 
form of a and uj mesons, which will subsequently decay primarily into pions with 
some photons also. The classical approximation is expected to be valid only when 
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Fig. 4: Calculated dependence of the center-of-mass radiated energy per nucleon upon impact pa- 
rameter for soft nucleon-nucleon collisions at five incident laboratory momenta. 

the amount of radiated energy in the center-of-mass system exceeds the mass of the 
hghtest meson, which is 550 MeV. As seen in fig. ^, this condition is well satisfied for 
impact parameters of physical interest at the three highest incident momenta, but not 
at the two lowest incident momenta. 

The qualitative behavior of these results can be understood in terms of the nature 
of the external forces. The repulsive vector force scales as the Lorentz factor 7 in 
both the longitudinal and transverse directions, whereas the attractive scalar force 
scales as 7^ in the longitudinal direction and as unity in the transverse direction. 
This implies that the vector force will dominate the transverse acceleration and the 
scalar force will dominate the longitudinal acceleration. For a given impact parameter 
the scattering angle and transverse momentum will be essentially proportional to the 
vector interaction strength (7^^, and the radiated energy will be essentially proportional 
to 7 times the scalar interaction strength Qs^. 



4. Future directions 

From nucleon-nucleon scattering experiments we know that several additional me- 
son fields are important and must be included for a quantitative description:ii 

• Isovector pseudoscalar (7r+, vr^, 7r°) 

• Isovector scalar {S~^, 6^, 5°) 

• Isovector vector (p+, p~, p°) 
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Fig. 5: Effective mass density for an extended nucleon. The intrinsic mass density is exponential with 
a root-mean-square radius of 0.862 fm. 

• Neutral pseudoscalar (rj) 

The next step in the systematic development of the theory should be the inclusion 
of these additional meson fields. Once this is done, the massive bremsstrahlung in 
our theory would include pions, deltas, rhos, and etas in addition to the sigmas and 
omegas that are produced in the current version. 

The effects of quantum uncertainty on the equations of motion should also be 
studied and included if they are important. This should be possible by use of tech- 
niques analogous to those used by Moniz and Sharp for nonrelativistic quantum 
electrodynamics.0 There appear in the classical nonrelativistic equations of motion 
for an extended electron terms of the form pi^XP p{^') where the oper- 

ator O is a function of r and r'. Moniz and SharpEJ have shown in nonrelativistic 
quantum electrodynamics that the effect of quantum mechanics on the equations of 
motion is to replace such terms by terms of the form p{r) Cgfr Pesij"') dh" cPr' and 

derivatives with respect to A of these terms, where A = l/(Mo) is the Compton wave- 
length associated with the particle's bare mass and the effective operator O^s is a 
function of r, r', and A^ Vr'^. As illustrated in fig. |^, the effective mass density for an 
extended nucleon oscillates around the intrinsic mass density as a function of radial 
distance from the origin.0 

In conclusion, we have shown that classical hadrodynamics requires minimal phys- 
ical input, leads to equations of motion that can be solved numerically without further 
approximation, and provides a suitable framework for describing particle production 
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in relativistic heavy-ion collisions through the mechanism of massive bremsstrahlung. 
This work was supported by the U. S. Department of Energy. 
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